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Abstract We present a study on low-T c superconductor- 
insulator-ferromagnet-superconductor (SIFS) Josephson 
junctions. SIFS junctions have gained considerable inter- 
est in recent years because they show a number of inter- 
esting properties for future classical and quantum com- 
puting devices. We optimized the fabrication process of 
these junctions to achieve a homogeneous current trans- 
port, ending up with high-quality samples. Depending on 
the thickness of the ferromagnetic layer and on temper- 
ature, the SIFS junctions arc in the ground state with a 
phase drop either or tt. By using a ferromagnetic layer 
with variable step-like thickness along the junction, we 
obtained a so-called 0-tt Josephson junction, in which 
and 7T ground states compete with each other. At a 
certain temperature the and tt parts of the junction 
are perfectly symmetric, i.e. the absolute critical current 
densities are equal. In this case the degenerate ground 
state corresponds to a vortex of supercurrent circulating 
clock- or counterclockwise and creating a magnetic flux 
which carries a fraction of the magnetic flux quantum 
<P . 



1 Introduction 

Superconductivity (S) and ferromagnetism (F) are two 
competing phenomena. On one hand a bulk supercon- 
ductor expels the magnetic field (Meissner effect). On 
the other hand the magnetic field for H > H C 2 destroys 
the superconductivity. This fact is due to the unequal 
symmetry in time: ferromagnetic order breaks the time- 
reversal symmetry, whereas conventional superconduc- 
tivity relies on the pairing of time-reversed states. It 
turns out that the combination of both, superconductor 
and fcrromagnet, leads to rich and interesting physics. 
One particular example - the phase oscillations of the su- 
perconducting Ginzburg-Landau order parameter inside 
the ferromagnet - will play a major role for the devices 
discussed in this work. 



The current-phase relation I s (4>) of a conventional 
SIS Josephson junction (J J) is given by I s {4>) = Ic sin(</>). 
(j> = 6\ — $2 is the phase difference of the macroscopic 



superconducting wave functions = 
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parameters of each electrode) across the junction, I c 
is the critical current. Usually I c is positive and the 
minimum of the Josephson energy U = Ej(l — cos</>), 
Ej = k^- is at 4> — 0- However, Bulaevskii et al. [1] cal- 
culated the supercurrent through a JJ with ferromag- 
netic impurities in the tunnel barrier and predicted a 
negative supercurrent, I c < 0. For —I c sm.(<p) = the 
solution <fi = is unstable and corresponds to the max- 
imum energy U = Ej(\ + cos(f>), while <fi = tt is sta- 
ble and corresponds to the ground state. Such JJs with 
(f> = 7r in ground state are called tt junctions, in contrast 
to conventional junctions with <j) — 0. In case of a tt 
Josephson junction the first Josephson relation is modi- 
fied to I s {4>) = — -^c sm (0) = Ic sin(4> + 7r). In experiment 
the measured critical current in a single junction is al- 
ways positive and is equal to |7 C |. It is not possible to 
distinguish JJs from tt JJs from the current-voltage 
characteristic (IVC) of a single junction. The particular 
I C (T) [2] and I c {d F ) [3] dependencies for SFS/SIFS type 
junction are used to determine the tt coupled state. For 
low-transparency SIFS junctions the I c {<If) dependence 
is given by 
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where £fij£f2 are the decay and oscillation lengths of 
critical current and d F ca,d is the dead magnetic layer 
thickness [4]. For \£,F2^ < c?f — d F cad < §£f27t the cou- 
pling in ground state of JJs is shifted by tt. 

In a second work Bulaevskii et al. [5] predicted the 
appearance of a spontaneous supercurrent at the bound- 
ary between a and a tt coupled long J J (LJJ). This 
supercurrent emerges in the absence of a driving bias 
current or an external field H, i.e. in the ground state. 
Depending on the length of the junction L the supercur- 
rent carries one half of the flux quantum, i.e. #o/2 (called 
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semifluxon) , or less. Fig. 1(a) depicts the cross section of 
a symmetric 0-tt long J J. The spontaneous supercurrent 
j s flows either clockwise or counterclockwise, creating 
the magnetic field of ±#o/2- The current density jumps 
from maximum positive to maximum negative value at 
the 0-7r phase boundary. A theoretical analysis based on 
the perturbed sine-Gordon equation is given in Ref. [6]. 
Below we will first discuss the properties of the sponta- 
neous supercurrent and, second, various systems having 
0-7r phase boundaries. 

Spontaneous supercurrent Kirtley et al. [7] calculated 
the free energy of 0-tt J Js for various lengths of the and 
tt parts as a function of the normalized length I = L/Xj 
and the degree of asymmetry A = \ j'^\L 7l /\j^.\LQ. where 
j^,jc are the critical current densities and Lq, L t are the 
lengths of and tt parts respectively, so that L = Lq+L w . 
The state of a symmetric 0-tt junction (A = 1) with 
spontaneous flux has lower energy than the states </> = 
or <f> = tt without flux. Symmetric 0-tt junctions have 
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Fig. 1 (a) Sketch of a 0-tt SIFS JJ with step-like thick- 
ness of F-layer and circulating supercurrent j s around 0-7r 
phase boundary. The junction length L > Aj, therefore the 
spontaneous flux (area below magnetic field) is equal to half 
of a flux quantum <3?o (semifluxon). (b)-(d) depicts the phase 
4>{x), magnetic field 4> x {x) and supercurrent j e (%) = jfj sin0 
of the 0-tv junction. 



always some self-generated spontaneous flux, although 
its amplitude vanishes for L — ► as <P w <Po£ 2 /8ir. For 
example, a symmetric 0-tt JJ of the total length L = Xj 
has a spontaneous magnetic flux <P « 0.04<?o and a sym- 
metric 0-7T JJ with L = 8Aj has a spontaneous flux of 
some 2 — 3% below Only in case of a infinitely long 

JJ we refer to the spontaneous flux as semifluxons, for 
shorter JJs it is named fractional vortex. 
The supercurrent or magnetic flux can be directly de- 
tected by measuring I C {H) [7], by scanning SQUID (su- 
perconducting quantum interference device) microscopy 
(in the LJJ limit, see [8, 9]) or by LTSEM (low temper- 
ature scanning electron microscopy) [10]. 



0-7T junctions technology 0-7T Joscphson junctions with 
a spontaneous flux in the ground state are realized with 
various technologies. The presence of fractional vortex 
has been demonstrated experimentally in e?-wave super- 
conductor based ramp zigzag junctions [9], in long Joseph- 
son 0-7T junctions fabricated using the conventional Nb/- 
Al-A^Oa/Nb technology with a pair of current injec- 
tors [11], in the so-called tricrystal grain-boundary LJJs 
[8, 12, 13] or in SFS/SIFS JJs [14, 15, 1G] with stepped 
ferromagnetic barrier as in Fig. 1. In the latter systems 
the Josephson phase in the ground state is set to or tt 
by choosing proper F-laycr thicknesses d\, e?2 for and 
tt parts, i.e. the amplitude of the critical current den- 
sities j° and can be controlled to some degree. The 
advantages of this system are that it can be prepared in 
a multilayer geometry (allowing topological flexibility) 
and it can be easily combined with the well-developed 
Nb/Al-Al 2 3 /Nb technology. 

The starting point for estimation of the ground state of a 
stepped J J is studying the IVCs and I C {H) for the planar 
reference and tt JJs. From this one can calculate im- 
portant parameters such as the critical current densities 
the Josephson penetration depths Aj, A} and the 
ratio of asymmetry A. For 0-7T junctions one needs and 
tt coupling in one junction, setting high demands on the 
fabrication process. The ideal 0-tt JJ would have equal 
be I = be I an d a _7r phase boundary in its center to 
have a symmetric situation. Furthermore the junctions 
should be undcrdamped (SIFS structure) since low dis- 
sipation is necessary to study dynamics and eventually 
macroscopic quantum effects. The junctions should have 
a high j c (and hence small Xj oc Vjc) to reach the LJJ 
limit and to keep high V c = I C R products, where V c is 
the characteristic voltage and R the normal state resis- 
tance. 

Previous experimental works on 0-7T JJs based on SFS 
technology [14, 15] gave no information about j° and ■ 
Hence, the Josephson penetration depth Aj could not be 
calculated for these samples and the ratio of asymmetry 
A was unknown. The first intentionally made symmetric 
0-7T tunnel JJ of SIFS type with a large V c was realized 
by the authors [16], making direct transport measure- 
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Fig. 2 (Color online) I C (H) of (a) SIFS (4 nm NiCu) and 
(b) SINFS (Cu 2 nm, NiCu 4.7 nm) stacks. Oxygen pressure 
is 0.45 mbar for SIFS and 0.015 mbar for SINFS type. 
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Fig. 3 (Color online) I c (cLf) and I C {T) (inset) dependences 
of SIFS junctions at 4.2 K. Note the difference of the slope 
of I C {T) for and tt coupled junction (inset). 



mcnts of I C (H) and calculation of the ground state with 
spontaneous flux feasible. 

Within this paper we review the physics of 0-7T cou- 
pled SIFS-type Josephson junctions and give an overview 
on our experimental results. Special focus is put on the 
fabrication of SIFS junctions having a planar or stepped- 
typcd ferromagnetic layer (NiCu), the determination of 
ground state (0 or tt for planar JJs) and asymmetry of 
critical currents (stepped JJs). Finally we give an esti- 
mation of the spontaneous magnetic flux in the ferro- 
magnetic 0-7T JJs. 

2 Fabrication 

The fabrication process for planar junctions is based 
on Nb/Al-Al 2 3 /NiCu/Nb stacks, deposited by dc mag- 
netron sputtering [17]. Thermally oxidized 4-inch Si wafer 
served as substrate. First of all, a 120 nm thick Nb 
bottom electrode and a 5 nm thick Al layer were de- 
posited. Second, the aluminium was oxidized for 30 min 
at room temperature in a separate chamber. Third, the 
ferromagnet (i.e. NigoCu4o alloy, Tc = 225 K) was de- 
posited. To have many structures with different thick- 
nesses in one fabrication run, we decided to deposit a 
wedge-shaped F-layer. For this the substrate and sputter 
target were shifted about half of the substrate diameter. 
This allowed the preparation of SIFS junctions with a 
gradient in F-layer thickness in order to minimize in- 
evitable run-to-run variations. The sputtering rates for 
NiCu along the gradient were determined by thickness 
measurements on reference samples using a Dektak pro- 
filer. At the end a 40 nm Nb cap layer was deposited. 
The tunnel junctions were patterned using a three level 
optical photolithographic mask procedure and Ar ion- 
beam milling [18] . The insulation between top and bot- 
tom electrode is done by a self- aligned growth of Nb 2 Os 
insulator by anodic oxidation of Nb after the ion-beam 
etching. The Nb20s exhibited a defect free insulation 
between the superconducting electrodes. 



Topological and electrical measurements, see Ref. [17], 
indicated that the direct deposition of NiCu on the tun- 
nel barrier (SIFS-stacks) led to an anomalous I C (H) de- 
pendence such as shown in Fig. 2(a), which is an indica- 
tion for an inhomogencous current transport. An addi- 
tional 2 nm thin Cu layer between the AI2O3 tunnel bar- 
rier and the ferromagnetic NiCu (SINFS-stacks) brought 
considerable benefits, as it ensured a homogeneous cur- 
rent transport, see Fig. 2(b). In this way a high number 
of functioning devices with j c spreads less than 2% was 
obtained. The variation of the F-layer thickness over a 
length of one junction diameter is less than 0.02 nm. For 
simplification we refer in the following to SIFS stacks, 
although the actual multilayer is SINFS-type. 
The patterning of stepped junctions was done after the 
complete deposition of the planar SIFS stack and before 
the definition of the junction mesa by argon-etching and 
Nb20s insulation. The detailed process is published in 
Ref. [19]. The junction was partly protected with pho- 
toresist to define the step location in the F-layer, fol- 
lowed by i) selective reactive etching of the Nb, ii) ion- 
etching of the NiCu by AoIf and iii) subsequent in situ 
deposition of Nb. To our knowledge, this was the first 
controlled patterning of 0-7T JJs based on a ferromag- 
netic interlayer. 

The planar 0, n reference junctions and the stepped 0-7T 
junctions were fabricated from a single trilaycr. 

3 SIFS junctions without step-like F-layer 

All investigated junctions had an area of 10 000/im 2 , but 
the length and width were different for different junc- 
tions. The length was comparable or shorter than the 
Josephson penetration depth A,/. We investigated the 
thickness dependence of the critical current I c (dp). To 
produce the AI2O3 barrier the Al layer was oxidized at 
0.015 mbar yielding j c as 4.0 kA/cm 2 for the reference 
superconductor-insulator-superconductor (SIS) JJs. Then 
SIFS stacks with wedge-like F-layer were fabricated in 
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another run. Taking the JJs of the same geometry (100 x 
100 /Jin 2 ), but situated at different places on the wafer 
(i.e. different <7f) we have measured the nonmonotonic 
I c {<If) dependence shown in Fig. 3. As a result the fit- 
ted parameters are £p± = 0.78 nm, £p2 = 1.35 nm and 
dje&d ~ 3 09 nm. The coupling changed from to tt at 
the crossover thickness dp"* = f £f2 + <i F oad = 5.21 nm 
[4]- 

The magnetic and spin-orbit scattering in the F-layer 
mixes the up and down spin states of electrons in the 
conduction bands. If the spin-flip scattering time t s is 
short ftr^r 1 3> fes^c, like in NiCu alloys, the temper- 
ature dependence of scattering provides the dominant 
mechanism for the I c (T) dependence [20] . The oscillation 
period £f2 becomes shorter for decreasing temperature, 
thus the whole I c (dp) dependence is squeezed to thinner 
F-layer thicknesses. Hence, the temperature dependence 
of the critical current 7 C (T) is an interplay between an in- 
creasing component due to an increasing gap and a mag- 
netic coupling dependent contribution which may de- or 
increase 7 C . The 7 C (T) relations for two JJs (one 0, one tt) 
arc shown in the inset of Fig. 3. At g?f = 5.11 nm the JJ 
is coupled, but one can relate the nearly constant 7 C be- 
low 3.5 K to the interplay between an increasing gap and 
a decreasing oscillation length £,F2(T). The dp = 5.87nm 
JJ is tt coupled and showed a linearly increasing I c with 
decreasing temperature. 

4 SIFS junctions with step-like F-layer 

Various structures on the wafer were placed within a 
narrow ribbon perpendicular to the gradient in the F- 
layer thickness and were replicated along this gradient. 
One ribbon contained reference JJs with the uniform 
F-layer thickness d\ (uniformly etched) and c?2 (non- 
etched) and a JJ with a step Adp in the F-layer thick- 
ness from di to di- The lengths Ld 1 and Ld 2 are both 
equal to 167 /mi. The lithographic accuracy is of the 
order of 1 yum. A set of structures with difference in 
dp between neighboring ribbons of 0.05 nm was ob- 
tained. Comparing the critical currents I c of non-etched 
JJs (dots), see Fig. 4 with the experimental I c (dp) data 
for the etched samples (stars) we estimate the etched- 
away F-layer thickness as Adp as 0.3 nm. The stars in 
Fig. 4 are shown already shifted by this amount. Now we 
choose the set of junctions which have the thickness di 
and critical current l c {d-i) < (tt junction) before etch- 
ing and have the thickness d\ = efe — AdF and critical 
current I c (di) as —I c {d2) (0 junction) after etching. One 
option is to choose the junction set denoted by closed cir- 
cles around the data points in Fig. 4, i.e. d\ = 5.05 nm 
and g?2 = 5.33 nm. 

The I-V characteristics and the magnetic field de- 
pendence of the critical current 7 C (77) was measured 
for all three junctions: J J with dp = d\, tt J J with 
dp = d-2 and 0-7T JJs with stepped F-layer (d\ and di 
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Fig. 4 (Color online) Critical current I c of the uniformly 
etched (star) and non-etched (dot) SIFS junctions versus the 
F-layer thickness before etching (If - The fit of the experimen- 
tal data for non-etched samples using Eq.(l) is shown by the 
continuous line. The JJs were oxidized at 0.015 mbar. 

in each half). The magnetic diffraction pattern 7 C (77) of 
the 0-7T JJ and the and tt reference JJs are plotted 
in Fig. 5. The magnetic field 77 was applied in-plane of 
the sample and parallel to the step in the F-layer. Due 
to a small net magnetization of the F-layers the I C (H) 
of references junctions were slightly shifted along the 77 
axis. Nevertheless, both had the same oscillation period 
Li a H cl as 36 f/T. At T as 4.0 K the 0-?r JJs was slightly 
asymmetric with 7° as 208 (xA and I™ as 171 fiA (data 
of reference JJs) . To achieve a more symmetric configu- 
ration, the bath temperature was reduced, because a de- 
crease in temperature should increase I£ = I c {d2) more 
than I® = I c (di), like for the and tt samples in the inset 
of Fig. 3. As a result, both 7°(T) and I£(T) were increas- 
ing when decreasing the temperature, but with different 
rates. At T as 2.65 K the critical currents 7° and 
became approximately equal, see Fig. 5. The magnetic 
field dependence of the planar reference junctions Ic(H) 
and Ic(H) look like perfect Fraunhofer patterns. One 
can see that the Ic(H) and 7^(77) measurements almost 
coincide, having the form of a symmetric Fraunhofer pat- 
tern with the critical currents 7° as 220 /jA, 7^ as 217 //A 
and the same oscillation period. The stepped 0-7T junc- 
tion had a magnetic field dependence 7°~'"'(77) with a 
clear minimum near zero field and almost no asymme- 
try The critical currents at the left and right maxima 
(146 /iA and 141 fiA) differ by less than 4 %, i.e. the 
0-7T junction is symmetric, and its ground state in ab- 
sence of a driving bias or magnetic field (7 = 77 = 0) 
can be calculated [16]. Our symmetric 0-7T LJJ had an 
normalized length of £ = 1.3, with a spontaneous flux in 
the ground state of 

±$ as <p f /8tt as 0.067 ■ <P , 
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Fig. 5 (Color online) I C (H) of 0-tt JJ (open triangles) with 
H applied parallel to short axis, overlayed with the non- 
etched (dot) and etched (stars) reference SIFS junction mea- 
surements. At T ~ 2.65 K the 0-tt JJ becomes symmetric. 
The junction dimensions are 330 x 30 /an 2 . 

being equal to 13% of <Po/2. A detailed calculation taking 
several deviations from the ideal short JJ model into 
account can be found elsewhere [21]. 

5 Summary 

The concept and realization of 0-tt junction based on 
SIFS stacks has been presented. The realization of tt cou- 
pling in SIFS junctions and the precise combination of 
and 7T coupled parts in a single junction has been shown. 
The coupling of the ferromagnetic Josephson tunnel junc- 
tions was investigated by means of transport measure- 
ments. The emergence of a spontaneous flux, which was 
calculated as 13% of half a flux quantum <Pq/2, was ob- 
served in the magnetic field dependence of the current- 
voltage characteristics of the 0-tt J J. 
As an outlook, the ferromagnetic 0-tt Josephson junc- 
tions allow to study the physics of fractional vortices 
with a good temperature control of the symmetry be- 
tween and tt parts. We note that symmetry is only 
needed for J J lengths L < Xj. For longer JJs the semi- 
fluxon appears even in rather asymmetric JJs, and T 
can be varied in a wide range affecting the semifluxon 
properties only weakly. The presented SIFS technology 
allows us to construct 0, tt and 0-tt JJs with comparable 
j° and j* in a single fabrication run. Such JJs may be 
used to construct classical and quantum devices such as 
oscillators, memory cells, tt flux qubits [22, 23] or semi- 
fluxon based qubits [24]. 
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